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Abstract 

We study a large population of globally coupled phase oscillators subject to common white 
s— / ' Gaussian noise and find analytically that the critical coupling strength between oscillators for 

r ' synchronization transition decreases with an increase in the intensity of common noise. Thus, 

common noise promotes the onset of synchronization. Our prediction is confirmed by numerical 
simulations of the phase oscillators as well as of limit-cycle oscillators. 
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Synchronization of an ensemble of periodic oscillators has attracted considerable attention 



IM 



because of its broad applications in many fields ranging from physics to engineering [iH^. In 
particular, synchronization plays an essential role in numerous biological functions, including 
the formation of pacemaker tissues of the heart and of the circadian master clock [Sj. 

Because real systems are inevitably subject to noise, it is important to understand the 
effect of noise on the synchronization of periodic oscillators. Some types of noise, including 
thermal noise or intrinsic noise in cells, act independently on individual components, which 
usually inhibits synchronization [l|, |6|. However, there are many situations where a single 
noise process, such as that originating from environmental fluctuations, acts on an entire 
system. Whether such common noise enhances or inhibits synchronization is actually un- 
clear. This issue is thought to be relevant to biological pacemaker tissues in that external 
noise could have a positive effect on synchronization. However, clarification of the outcome 
of fluctuating input is necessary in cases, such as that of deep brain stimulation for Parkin- 
son's disease [7|, in which global external stimulation is used to destroy synchronization of 
dynamic components. 

The effect of common noise on uncoupled oscillators or coupled-oscillator networks with 
small sizes has been extensively studied for both periodic and chaotic oscillators, and rigorous 
theoretical frameworks have been proposed 8|. In contrast, for a large population of coupled 
oscillators, despite a numerous body of numerical and experimental evidence [9| , theoretical 
treatment is still an open and challenging problem. 

In this letter, we investigate a large population of nonidentical phase oscillators that 
are globally coupled and subject to common Gaussian white noise to clarify the effect of 
common noise on coupled oscillators. Utilizing the anzatz recently proposed by Ott and 
Antonsen [lO(], we analytically show that the addition of common noise leads to a decrease in 
the critical coupling strength for synchronization transition. Our prediction is corroborated 
by direct numerical simulations of the model. We also numerically confirm that globally 
coupled limit-cycle oscillators show the same dependence on common noise. The employed 
phase model approximates many realistic systems with weak coupling and weak forcing. 
Thus, our results suggest that weak common noise generally promotes synchronization of 
oscillators with week and global coupling. 

Consider globally coupled phase oscillators, known as the Sakaguchi-Kuramoto model 



ll| . subject to a common external force 



df) K ^ 

d^ = ^^ + iV S '''^^^' - ^^ + /3) + V{t) sin Oi, (1) 

i=i 

where 6*4 and Ui are the phase and the natural frequency, respectively, of the i-th oscillator, 
f^ > is the coupling strength, /3 is a parameter of the coupling function (— vr/2 < [3 < t^ 12), 
and p{t) is a common external force. We assume that the natural frequency distribution is 
given by a Lorentzian function /froq(w) = \ i^^_^ \2,i - We will further assume white Gaussian 
noise for p{t) [12], but first we treat p{t) as a general time- dependent function for a while. 



Note that 
forcing [l|, y. 



. dl]) approximates various realistic oscillators with weak coupling and weak 
13[. Note also that the common external force is multiplied by a function of 
the phase, sin6'j, which is called a phase sensitivity function. The phase sensitivity function 
naturally appears in the phase description of limit-cycle oscillators [l|, |2|, [l3|. We will later 



demonstrate these facts using a limit-cycle-oscillator model that generally appears near a 
Hopf bifurcation. 

We examine the synchronization transition in the large- A^ limit. For a better presentation, 
we set /3 = (corresponding to the Kuramoto model [l]). The extension to nonzero /3 is 
straightforward; we will only show a final result for nonzero /3 in the present paper. In the 
limit A^ — 7- 00, Eq. ([1]) becomes 

where f{u,6,t) is the distribution function for the phases of the oscillators with natural 
frequency u, r = J_^du J^ ^d6 /e'^ is the Kuramoto order parameter, and * represents the 
complex conjugate. For p{t) = 0, the synchronization transition (the so-called Kuramoto 
transition) occurs at K = K^ = 2, above which \r\ is nonvanishing [l|. 

To investigate the transition in Eq. ([2]), we first derive a dynamical equation for the order 
parameter r. For this, we employ the Ott-Antonsen ansatz for the phase distribution 



/ = ^^|l + E[(-T+(«*e-)"]|, (3) 

where a{u,t) is a certain function [lO|. Substituting Eq. ([3]) into Eq. (|2]), we obtain a 
dynamical equation for a 

^ + y (m^ - n + ica + |(1 - a') = 0. (4) 



Note that r(t) = a*{ijJo — i, t) because r = J^ do; J^ '" d^/e'^ = /^ duff^^qC^* and ffj-cq 
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TV. Thus, by setting to = loq — iin Eq. (jl]), we obtain 



^=(-l + y + -ojr--|r|V--(l-r^). (5) 

By letting r = y/Ae^^'^°*^'^\ we further obtain 

dA 

— = h{A) + gA{A, uot + (f))p, (6) 

— = g4,{A,uJot + (f))p, (7) 

where h{A) = {K -2)A-KA'^, gA{A,u}ot + (l)) = -y/A{l- A) cos{ujot + <!)), and g^{A,ujot + 

0) = ^sin(u;ot + 0). 

Now we assume that p{t) is white Gaussian noise with {p{t)) = and {p{t)p{s)) = 
2DS{t — s), and interpret Eq. ([7]) as a Stratonovich differential equation. Then we obtain 
the Fokker-Planck equation for the probability distribution q{A,(j),t), given by 

dq 9 If f dgA ^ dgA\\ \ ^ ( f dg^ dgA 

Because h, gA, (?</,, and q are 27r-periodic functions, integrating of both sides of Eq. ([8]) over 
(J) from to 27r yields 

where Q(A, t) = j^^ q dcp. 

At this stage, we additionally assume that K and D are sufficiently small compared to 
a typical natural frequency ojq. It is natural to assume this because this is the condition 
under which Eq. ([T]) approximates coupled limit-cycle oscillators. Under this assumption, Q 
evolves sufficiently slowly compared to a typical oscillation time scale, i.e., 2t\: /ojq. Thus, to 
a good approximation, the right-hand side of Eq. ([9]) can be time-averaged over the duration 
of 27r/(X'o, leading to 



dt M i V 2 V 2 J J ^ } dA^ [2 

Letting dQ/dt = 0, we obtain the stationary distribution Qoo{A) as 

' 2 r 2A 
Qoc(A) = Cexp -}--—^-(K + D)\og{l-A) 



(11) 
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FIG. 1. (color online) Numerical results for the phase model given by Eq. ([T]). Crosses (black) and 
open circles (orange) represent numerical data for D = and D = 0.02, respectively, (a) Snapshot 
of the phase distribution for K = 1.99. (b) Distribution of A for (b-1) K = 1.96, (b-2) K = 1.99, 
and (b-3) K = 2.1. Lines on the points are fitting curves. Histograms and curves are normalized 
for the maximum of curves to be 1. Point-dashed line (blue) and dashed line (green) represent the 
numerically identified ^max for D = and D = 0.02, respectively. 



where C = 1/ J^ dAexp [-^ {~T^ ~ (-^ + ^) los(l ~ ^)}] • -'■^ stochastic systems, the max- 
iinuin of the probabiUty distribution function is often adopted as the order parameter char- 
acterizing a transition [1^]. From Eq. ( ITTl) . it follows that Qoo(^) assumes its maximum 
at 
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Thus we find that the critical coupling strength at which Amax becomes nonvanishing is 
Kc = 2 — D; the common noise decreases the critical coupling strength by D as compared 
to that in the original Kuramoto transition. 

For nonzero (3, one can show that K in Eq. ( 1T2|) is replaced by K cos/3. Thus, the critical 
condition is given by K^ 



2-D 

cos/3 ■ 



We confirmed our prediction by numerical simulation of Eq. ([T]) with A^ = 10000 and 



P = 0. The Lorentzian distribution for the natural frequency was given by ISj 

a;, = a;o + tan|*^-(iV + l)^} (1 < z < iV). (13) 

We set Wo = 100 to ensure that K and D are much smaller than uq- We employed random 
initial conditions and numerical data were obtained from t = 10000 to t = 60000. As shown 
in Fig. [T](a), the phase distribution did not cluster for K = 1.99 and D = [K + D < 2). In 
contrast, a cluster of oscillators was observed for K = 1.99 and D = 0.02 {K + D > 2). To 
estimate Amax from the numerical data, the logarithm of the histogram of A around the peak 
was fitted to the logarithm of Eq. ([IT]), i.e., a + | (- j^ - (c + b) log(l - A)) with fitting 
parameters a, b, and c. The obtained data were well fitted [Fig. [1] (b)]. The numerically 
identified values of Amax = ^~hTc^ were plotted in Fig. HI which shows excellent agreement 
with the theoretical prediction of Eq. ( TT2|) . In our preliminary numerical simulations, we also 
confirmed that a similar transition behavior occurs in the case of the Gaussian distribution 
for the natural frequency (data not shown). 

We also observed the distribution of the averaged frequencies cjf™ defined as the long-time 
average of 6i. Numerical results are shown in Fig. [3l Without noise, the distribution had 
a delta-function peak at Uq, whereas for D ^ 0, this peak disappeared and the distribution 
was continuous. This qualitative difference can be explained as follows. Using the Kuramoto 
order parameter r, Eq. ([1]) can be written as 9i = Ui + K\r\ sin(u;o'^ + (p — ^i) +p{t) sin 9i. For 
D = 0, \r\ and are time-independent after transient [l]. Then, oscillators with \ui — ujqI < 
\r\ are phase-locked to the mean field, so that their actual frequencies are exactly the same 
as that of the mean field, which is uq. However, for D ^ 0, \r\ fluctuates with time and 
becomes vanishingly small with a finite probability [see Eq. ( ITTl) and Fig.[T](b)]. This implies 
that any oscillator except that with Ui = uq cannot be phase-locked to the mean field for 
an infinitely long time. Therefore, oscillators with Ui > Uq (wj < Uq) tend to have a larger 
(smaller) averaged frequency than that for D = 0, so that the delta-function peak vanishes. 

Finally, we demonstrate the validity of our prediction in limit-cycle oscillators. We intro- 
duce the following model 

r\W K ^ 

^ = (1 + iuj,)W, - \W,\^ W^^ + ^ E(^^- - ^^) + v^P(^)' (14) 

where Wi is the complex state variable of the z-th oscillator, e is a small parameter to denote 
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FIG. 2. (color online) Amax for the phase model as a function of K. Crosses (black), open circles 
(orange) and filled circles (red) indicate the numerically identified A^yax for D = 0, D = 0.02, and 
D = 0.04, respectively. Error bars represent the variance of j4max for 10 trials with different initial 
conditions and different noise processes. Point-dashed line (black), dashed line (orange), and a 
solid line (red) represent Eq. (|12p for D = 0, D = 0.02, and D = 0.04, respectively. 
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FIG. 3. (color online) Distribution of the long-time averaged frequencies of the phase oscillators 
for K = 2.02. Crosses (black) and open circles (orange) with connecting lines are the numerical 
results for D = and D = 0.02, respectively. 



that the coupling strength and the noise strength smaller than both the relaxation rate of 
the amplitude dynamics and the natural frequencies of oscillators, and p{t) is a common 
white Gaussian noise with strength D. Each individual oscillator is called a Stuart-Landau 
oscUato., which ge„e*a..y appea. when .he sys.e. is nea. a Hopf b,fu.ca«o„ Q. E,. O 
is approximated by Eq. ([1]) with (3 = for small e [1|, so similar behavior is expected. 
We numerically simulated Eq. dH]) with A^ = 1000. We defined A as lE^^ie'^V^P 
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FIG. 4. (color online) Numerical results for the limit-cycle oscillators given by Eq. (|14p . Legends 
for (a) and (b) are the same as those in Figs. [2] and [H respectively, (a) Order parameter ^max as 
a function of K. Lines represent Eq. (J12p . (b) Distribution of the long-time averaged frequencies 
for K = 2.02. We used the data from t = 5 x 10^ to i = 10 x 10^. N = 1000, e = 0.01, and 
0.1 + etan{ifj-{N + l)^}. 
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with 9j = arg Wj and estimated Amax in the same manner as for the phase oscillators. The 
numerically determined A^ax values are shown in Fig. H] (a), which agrees reasonably well 
with the prediction of Eq. ( 1121) . We also observed that the distribution of uf"^ was continuous 
ioT D^O [Fig. 11(b)]. 

To conclude, we have studied the Sakaguchi-Kuramoto model subject to common 
noise and analytically showed that the critical coupling strength for the synchronization- 
desynchronization transition decreases with an increase in the strength of the common noise. 
The prediction has been numerically corroborated. We have also found that the distribution 
of the averaged frequencies is continuous when common noise is present. Our results suggest 
that weak common noise generally promotes synchronization of weakly coupled oscillators. 
It would be interesting to experimentally investigate the effect of common noise on coupled 
biological and chemical oscillators. 

We thank Dan Tanaka for motivating us to study this topic. We also thank Hayato Chiba 
and Hiroya Nakao for helpful discussions. 
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